Vortex and translational currents due to broken time-space symmetries 
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We consider the classical dynamics of a particle in a d = 2, 3-dimcnsional space-periodic potential 
under the influence of time-periodic external fields with zero mean. We perform a general time- 
space symmetry analysis and identify conditions, when the particle will generate a nonzero averaged 
translational and vortex currents. We perform computational studies of the equations of motion and 
of corresponding Fokker-Planck equations, which confirm the symmetry predictions. We address the 
experimentally important issue of current control. Cold atoms in optical potentials and magnetic 
traps are among possible candidates to observe these findings experimentally. 
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The idea of directed motion under the action of an ex- 
ternal fluctuating field of zero mean goes back to Smolu- 
chowski and Feynman It has been intensively stud- 
ied in the past decades again It is beheved to be 
connected with the functioning of molecular motors, and 
can be applied to transport phenomena which range from 
mechanical engines to an electron gas (see d, Q and ref- 
erences therein). 

The separation of the fluctuating fields into an uncorre- 
lated white noise term and a time-periodic field was used 
to perform a symmetry analysis of the most simple case - 
a point-likeparticle moving in a one-dimensional periodic 
potential [5|. It allowed to systematically choose space 
and time dependencies of potentials and ac fields such, 
that a nonzero dc current is generated. Various studies 
of the dynamical mechanisms of rectification have been 
reported (e.g. Q). Among many experimental reports, 
we mention the successful testing of the above symme- 
try analysis using cold atoms in one-dimensional optical 
potentials [2]. By use of more laser beams, experimen- 
talists can already fabricate two- and three-dimensional 
optical potentials, with different symmetries and shapes 
[8| , with the aim of even more controlled stirring of cold 
atoms in these setups. 

A particle which is moving in a d = 2, 3-dimensional 
periodic potential may contribute to a directed current 
along a certain direction. At the same time, the par- 
ticle can perform vortex motion (which is not possible 
in a one-dimensional setting) generating a nonzero av- 
erage of the angular momentum. Directed translational 
currents are supported by unbounded trajectories while 
vortex currents may be localized in a finite volume. The 
question is then, how can we control a type of the di- 
rected motion? To answer this question, one has to find 
conditions for an appearance of either purely vertex or 
purely translational currents. This is the main goal of 
the present work. 



In this letter we perform a symmetry analysis of par- 
ticle motion in a d = 2, 3-dimensional periodic potential, 
under the influence of external ac fields. We identify the 
symmetries which ensure that either directed currents, 
or vortex currents, are strictly zero. Breaking these sym- 
metries one by one allows to control the particle motion, 
to generate either directed, or vortex, currents, or both 
simultaneously. 

We consider the dynamics of a classical particle (e.g. 
an atom of a cold dilute atom gas, loaded onto a proper 
optical lattice) exposed to an external potential field: 

/^r + 7r = g(r,t)+^(i) , g(r, t) = - V;7(r, t) . (1) 

Here r — {x,y,z} is the coordinate vector of the par- 
ticle, the parameter 7 > characterizes the dissipation 
strength, and /i > defines the strength of the inertial 
term [9]. The force g(r,t) = {^^(r, i)}, a = x,y,z, is 
time- and space-periodic: 



g(r,t) =g(r,t + r) =g(r + L,i). 
The absence of a dc bias implies 



(2) 



L.T 



;(r, t) dt dxdydz — (3) 



where the spatial integration extends over one unit cell. 

The fluctuating force is modeled by a (5-correlated 
Gaussian white noise, $,{€) = {£,x,£.y,£,x}, {^a{t)^f3{t')) = 
2jDS{t — t')6aj3 {a,f} — x,y,z). Here D is the noise 
strength. The statistical description of the syste m (fll) is 
provided by the Fokker-Planck equations (FPE) : 

P(r,v,t) = {-Vr-v + Vv-[7v-g(r,t)] + 

-t- 7i?Av}P(r, V, t), M = 1, V = f , (4) 
7P(r,t) =- [Vr-g(r,<) + Z?Ar]P(r,t), Ai = (5) 

Each of the linear equations ([4][5l) has a unique attractor 
solution, P which is space and time periodic 11 1 . 



Directed transport. Let us consider the dc component 
of the directed current = v = r in terms of the 
attractor P: 

J= (vP(r,v,i))T,L, M = 1 , (6) 

J = i(g(r,t)-P(r,t)}T,L, M = 0. (7) 

7 

The strategy is now to identify symmetry operations 
which invert the sign of j, and, at the same time, leave 
Eq. ([T]) invariant. If such symmetries exist, the dc cur- 
rent J will strictly vanish. Sign changes of the current can 
be obtained by either inverting the spatial coordinates, 
or time (simultaneously allowing for shifts in the other 
variables). Below we list all operations together with the 
requirements the force g and the control parameters have 
to fulfill: 



Si : r 
^2 : r 

^2(g)^g (if7-0); 52(g) 



. -r + r', t ^t + r; 
r + A, t ^ -t + t'; 



g 



(if /i = 0). 



(8) 
(9) 



Here t' and r' depend on the particular shape of g(r,t). 
The system must be invariant under a spatial translation 
by the vector 2A in space and 2t in time, respectively. 
The vector A is therefore given by A = naJ-ia/^, 
Ua = 0,1, while r = 0,T/2. By a proper choice of g 
all relevant symmetries can be broken, and one can then 
expect the appearance of a non-zero dc current J [12^ . 

To be more precise, we consider the case of a particle 
moving in a two-dimensional periodic potential and being 
driven by an external ac field: g(r,<) = — VF(r)+E(i) = 
i{r) + 'E{t). The symmetry Si holds if the potential force 
is anti-symmetric^ f(~r + r') = — f(r), and the driving 
function is shift-symmetric, E(i -|- T/2) = — E(i). The 
symmetry 5*2 holds at the Hamiltonian limit, 7 = 0, if the 
driving force is symmetric, E(t-\-t') = E(— t). Finally, the 
symmetry ^2 holds at the overdamped limit, /i = , if the 
potential force is shift-symmetric, f (r -I- A) = — f (r) and 
the driving force is anti-symmetric, E(t + 1') = — E(— t). 

In order to break the above symmetries, we choose 

V{r) = V{x, y) = cos(a;)[l -I- cos(22/)] , (10) 
E^^yit) = Ei]l sint + Ei^l sm{2t + 9) . (11) 

The potential pH)) is shift-symmetric, A — {±7r,0}. The 
symmetry Si is broken since E is not shift-symmetric. 
Therefore in general we expect J ^ 0. 

In Fig[T]we show the computational evaluation of equa- 
tions (|4l5p |13| . We confirm the presence of a nonzero dc 
current. Applying operations Si and 9 ^ 9 + t: we con- 
clude J{9+Tr) = —J{9), which allows for an easy inversion 
of the current direction, as also confirmed by the data in 
FigUJa). In the overdamped hmit fj, = 0, S2 is restored 
for 9 = 0,±7r, and therefore J{—9) = —J{9) (thick lines 
in FigUJa)). Upon approaching the Hamiltonian limit. 
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FIG. 1: (color online) (a) Dependence of the current compo- 
nents, (solid line) and Jy (dashed line), on 9 for (fT|). (|10[l - 
(HH), with D = 1, Ei'-'> = -E^^ = 2, = -^^f' = 2.5. 

Data are for the overdamped (/i = 0,7 = 1, thick lines) 
and underdamped (/i = 1,7 = 0.1, thin lines) cases, respec- 
tively; (b) The time evolution of the mean particle position, 
f{t) = J rP(r, V, t)drdv, for 6 = 0. The trajectories are su- 
perimposed on the contour plot of the potential (I10|) . Curve 



(i) corresponds to iji^' = 3, E., 



r(2) 



E. 



(1) 



0, Ei'^ 



3.5 



and curve (ii) - to the parameters of panel (a). The other 
parameters are: D = = 1,7 = 0.1; (c) The phase lag si"' 
as a function of the dissipation strength 7. 



7 — > 0, the points where J = shift from 6* = 0, tt to 
9 = ±7r/2 where the symmetry S2 is restored (thin lines 
in Fig[TJa)). In the underdamped regime, the dc-current 
can be approximated as Ja oc Ja"^ sm[9 — 9a\'j)], a — 
x,y. The phase lag is equal to 6*^"^ — tt/2 and 9^Jj = 
in the Hamiltonian and overdamped limits, respectively 
(Figdc)) [M]. 

Even more control over the current direction is possi- 
ble, by imposing the symmetry conditions ([HllHl) on each 
component ga{y^, t) independently. For PH)) and pTjl with 



e: 



(2) 



(1) 



0, 9 



the symmetry transformation 
Sc '■ X —>■ —X, y y, t t-\-TT implies that the cur- 
rent along the x-direction is absent, — 0, and directed 
transport is happening along the y-axis [see Fig. lb, curve 
(i)] . We may conclude, that the symmetry analysis turns 
out to be a powerful tool of predicting and controlling 
directed currents of particles which move in d = 2,3- 
dimensional potentials under the influence of external ac 
fields. Note that dynamical mechanisms of current recti- 
fication of a two-dimensional deterministic tilting ratchet 



were discussed in Ref. [16 1. 

Vorticity. At variance to the one-dimensional case, 
particles in two and three dimensions can perform vor- 
tex motion, thereby generating ring currents, or nonzero 
angular momentum. First of all we note that the particle 
dynamics is not confined to one spatial unit cell of the 
periodic potential U{r,t). Even in the case when a di- 
rected current is zero due to the above symmetries, J = 0, 
the particle can perform unbiased difi^usion in coordinate 
space. In order to distinguish between directed transport 
and spatial diffusion on one side, and rotational currents 
on the other side, we use the angular velocity [l3| 



(12) 



as a measure for the particle rotation, where (...}t — 
limt^oo J Jq ■■■dt' . Q,{t) is invariant under translations in 
space and time. It describes the speed of rotation with 
which the velocity vector f (the tangential vector to the 
trajectory r(i)) encompasses the origin. 

Using the above strategy, we search for symmetry op- 
erations which leave the equations of motion invariant, 
but do change the sign of the angular velocity. If such 
symmetries exist, rotational currents strictly vanish on 
average. The sign of can be inverted by either (i) 
time inversion t ~t together with an optional space 
inversion r ±r, or (ii) the permutation of any two 
variables, e.g. Vxy ■ {x,y,z} {y,x,z}. That leads to 
the following possible symmetry transformations: 

R, -.r^Vr + r', t ^ t + t ; R,{g) ^ g , (13) 

i?2(g) (if7 = 0); i?2(g)^-g(if/i = 0). (14) 

Here V stands for any of the following operations: Vxy, 
Vyz or Vzx and t' and r' again depend on the particular 
shape of g(r, t). 

To be concrete, we will again consider a particle mov- 
ing in a two-dimensional periodic potential and being 
driven by an external ac field: g(r,i) = — Vy(r)+E(t) = 
f (r) -I- E(t). For d— 2 there is an additional transforma- 
tion due to a mirror reflection at any axis, S^, : {x, y} 
{a;, -y} or ty : {x, y} {-x, y}, 



g 



(15) 



Symmetry Ri can be satisfied for the Hamiltonian, un- 
derdamped and overdamped cases if Vi{Vr + r') = f(r) 
and 'PEit + t') = E(i). The symmetry R2 apply for the 
same cases as their counterparts S2 if there is no space 
inversion. In the presence of space inversion they can be 
satisfied both in the Hamiltonian [if f(— r) = — f(r) and 
E(t + t') = -E(-t)] and overdamped [if f (-r) = f (r) 
and E(t + t') = E(— f)] hmits. The symmetry R3 is 
relevant if can be applied: fx(x,—y) — fx{x,~y), 
fy{x,-y) = -fy{x,-y), Ex{t + T/2) = Exit) and 




FIG. 2: (color online) (a) Dependence Jni9), Eq.^, for 
©, Gil-dni, with fi ^ 1, D ^ 0.5, eP = 0.4, = 0.8 
and 7 = 0.2 (solid line), 7 = 0.05 (dashed line), and 7 = 2 
(dashed-dotted line). Insets: the trajectory (left insert) and 
the corresponding attractor solution, r{t), (right inset) for the 
case 7 — 0.2 and 9 — 7r/2. We have used A'' = 10^ independent 
stochastic realizations to perform the noise averaging; (b) The 
phase lag 60 as a function of the dissipation strength 7. 



Ey{t + T/2) = —Ey{t). Similar conditions can be found 
for Sj,. 

We performed numerical integrations of the equation 
of motion ([1]) with the following potential and driving 
force: 

V{x,y) — [—3 (cos.T -I- cosy) + cosxcosy]/2, (16) 
Ex{t) = Ei^^ cost , Eyit) = cos(t + 6) . (17) 

Averaging was performed over = 10^ different stochas- 
tic realizations 18]. Fig[2] shows the dependence of the 
rotational current ([12]) on the relative phase 9. The sys- 
tem is invariant under the transformation 5*1 ([8]), there- 
fore the directed current J = 0. However, for the under- 
damped case, 7^0, all the relevant symmetries (1131115^ 
are violated, and the resulting rotational current ([T^ 
is nonzero, and depends on the phase 6 (Fig. 2a). Note 
that symmetry R2 is restored when 6 — 0, ±7r, thus the 
current disappears in the Hamiltonian and overdamped 
limits for these values of the phase. The left upper inset 
in Fig[2| shows the actual trajectory of a given realiza- 
tion, confirming that the particle is acquiring an average 
nonzero angular momentum, while not leaving a small fi- 
nite volume due to slow diffusion and absence of directed 
currents. 

The exact overdamped limit, fj, — 0, is singular for the 
definition since the velocity of a particle, r(t), is a 
nowhere differentiable function. The overdamped limit 
can be approached by increasing 7 at a fixed /i = 1. Al- 
ternatively, one may remove the restriction on /i allowing 
for an infinitesimal value < /x <IC 1 at a fixed dissipation 
strength 7=1. Both parameter choices equally regular- 
ize . Numerical simulations for the former way of reg- 
ularization show that if 7//X > 5 then the rotational cur- 
rent completely refiects the symmetries corresponding to 



4 



the overdamped case (see dependence ^0(7) in FigllUb)). 

Let us discuss the relation of our results to the case 
of multidimensional stochastic tilting ratchets under the 
influence of a colored noise studied previously [l^ . Since 
equivalent (in a statistical sense) stochastic processes, 
^{t), have been used as driving forces, the symmetry Ri 
(fT5)l can be violated only by an asymmetric potential. 
But all potentials considered in Refs. [3| are invariant 
under the permutation transformation V. As a conse- 
quence, vortex structures for a local velocity field pre- 
sented in Refs. are completely symmetric (clockwise 
vortices are mapped into counterclockwise ones by V) 
and, therefore, the average rotation for any trajectory 
equals zero. 

The phase space dimension is five for d = 2 and seven 
for d = 3. Therefore, in the Hamiltonian limit (7 = 0), 
Arnold diffusion [20*1 takes place. The particle dynamics 
is no longer confined within chaotic layers of finite width. 
That leads to unbounded, possibly extremely slow, diffu- 
sion in the momentum subspace via a stochastic web [20| . 
Therefore a direct numerical integration of the equations 
of motion may lead to incorrect conclusions. 

There is rich variety of physical systems, where multi- 
dimensional ratchets can be observed: cold atoms in two- 
and three-dimensional potentials (optical guiding) 21j, 
colloidal particles on magnetic bubble lattices [22|, fer- 
rofluids [23|, and vortices in superconducting films with 



pinning sites 2J|. Our methods of directed current con- 
trol can be used for an enhancement of the particle sep- 
aration in the laser beams of complex geometry pBj . 

To conclude, we formulated space-time symmetries for 
the absence of both directed currents, and rotational cur- 
rents, for particles moving in spatially periodic poten- 
tials, under the influence of external ac fields. Proper 
choices of these potentials and fields allow to break the 
above symmetries, and therefore to generate and control 
directed and rotational vortex currents in an indepen- 
dent way. Numerical studies supplement the symmetry 
analysis and confirm the conclusions. 
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